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Abstract. 

Diffraction in time of matter waves incident on a shutter which is removed at time 
< = is studied in the presence of a linear potential. The solution is also discussed in 
phase space in terms of the Wigner function. An alternative configuration relevant to 
■ current experiments where particles are released from a hard wall trap is also analyzed 

' for single-particle states and for a Tonks-Girardeau gas. 

m 
o 

O ; PACS numbers: 03.75.-b, 03.75.Be, 0.3.75.Kk 

^ 1. Introduction 

s 

^ I Diffraction in time (DIT) was studied first by Moshinsky [1, 2], and a fiurry of botli 
^1 experimental and theoretical work investigating similar transient effects and setups have 
been carried out ever since. (For a brief, recent review see [3].) Quantum temporal 
^ ' oscillations of matter waves released from a shutter or confinement region constitute the 
' hallmark of the effect. 

Remarkably, all the theoretical studies have dealt with the quantum dynamics in 
free space or potentials with finite support [4, 5, 6]. However, the presence of a linear 
potential is includible for some experimental setups. This is actually the case for the 
first experimental observation of DIT, which made explicit use of gravity on a system 
of ultra cold atoms [7, 8] . Gravitational effects are also crucial in cold atom fountains 
as those used for frequency standards [9] . In addition, a linear potential arises from the 
interaction of a charged particle with an electric field, but it is within the field of atom 
optics that linear potentials, generally created by a magnetic field, have received a great 
deal of attention for ultracold atom manipulation, as a plugging potential for loading 
schemes, or aimed to cloud focusing, cooling, and other key operations [10, 11]. 

One of the natural configurations to deal with such situations in the stationary 
regime, namely that of a point-like and Gaussian sources, has already been discussed 
in excellent agreement with photo-detachment and atom-laser experiment [12]. In this 
work we generalize the Moshinsky shutter problem, i.e., a cut-off plane wave released at 
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time t = 0, by turning on a linear potential at that instant. The exact solution is derived 
in section 3 together with its Wigner function. A related configuration, namely, that of 
particles initially trapped in a rectangular box is analyzed in section 4. The experimental 
build up of all-optical hard-wall traps [13] has excited a great deal of attention for such 
geometries both at single [14, 15] and many-particle level [16, 17, 18, 19]. In particular, 
experiments are in view to study the quantum dynamics in the Tonks-Girardeau regime 
with Fock states of low particle number, N [20]. The experimental observation of 
DIT with a Bose- Einstein condensate [21] is in excellent agreement with the dynamics 
entailed in the time-dependent Schrodinger equation, showing that the condensate 
evolves as a free non-interacting particle gas after few millisecond of expansion. Indeed, 
the importance of the non-linearity during free expansion has been studied in full detail 
for the case of a Tonks-Girardeau gas in [19] and shown to dominate in a time scale 
shorter than tjsf = mL"^ /{2N7rh), L being the size of the initial trap and m the mass of 
the particle. The expansion of a Tonks-Girardeau gas will be analyzed in section 5. 



2. Diffraction in time in free space 

The paradigmatic setup for DIT corresponds to a beam of particles incident on a shutter 
represented by an infinite potential barrier. The relevance on the diffraction pattern of 
the reffectivity of the barrier [1] as well as different types of beam have been thoroughly 
investigated in [3]. For simplicity, we will assume here a totally absorbing barrier so 
that the initial state is given by a cut-off plane wave, 

^{x,t = 0) = e'P^/^Q{-x). (1) 

The solution to the time dependent Schrodinger equation [1, 2] reads 

e 2th 

ip{x,t) = M{x,p/h,ht/m) := — - — w{—z), (2) 



with 

l + i I t 



(p-^), (3) 



2 V mh V t 
and the so called Faddeyeva function [22, 23] w is defined as 



w(z) :— e ^ eric{—iz) — — du , (4) 

in Jy_ u — z 

where r_ is a contour in the complex z-plane which goes from — oo to oo passing below 
the pole. After [24, 1], M{x, k, t) has been named the Moshinsky function. The hallmark 
feature of this phenomenon is that, at variance with the classical solution, the probability 
density presents genuine quantum oscillations in time and space. Notice that a point of 
constant probability P propagates as 



, , pt nht , . 

xp{t) = \ Mo,p, (5) 

m \ m 

Uq^p being a given constant. This result will be compared with the case where the 
evolution takes place in the presence of a linear potential. 
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3. Diffraction in time with a linear potential 

We next study the DIT in the presence of a hnear potentiaL Consider a beam of particles 
incident from the left as before on a totally absorbing shutter located at position a; = 0, 
see (1). Suddenly, the shutter is removed at time equal zero, and a linear potential is 
switched on. Such model can simulate the quantum dynamics of a confined charged 
particle, with a homogeneous electric field turned on at the time the infinite potential 
wall at the origin is removed. Also, it can be applied to a particle in a gravitational 
field; in that case the cut-off plane wave may be used as a basis element to describe the 
actual initial state. 

The evolved wave function obeys the integral equation 

/oo 
dx'K{x,t\x',0)tjj{x',0), (6) 
-oo 

where K{x,t\x',t') is the propagator for a Hamiltonian with a general linear potential 
of the form 

H = ^^f-, (7) 
with / constant. It is indeed well known [25], 

K{x,t\x',t') = J—— e 2n 24a ^ (g) 



27rih(t - 1') 

and related to the action of the classical path. 
Inserting (8) into (6) we find 



m 







It is convenient to introduce the auxiliary variables 



^(^x,t) = \l——e-'^-'^ I dx'e'^nr^+'^-^" (9) 



il){x,t) = -^e ^'a^ I due'"^ , (11) 

oo 



to write (9) as 

which can be expressed in terms of the Presnel sine and cosine integrals, 

/"() 2 1 _|_ 7 /""O 2 1-1-7 

due'^^^+ due'^ ^^ + C{uo)+iS{uo). (12) 
oo 2 7q 2 

Prom the last equation it is possible to map the probability density to the Cornu spiral 
[1]. Additionally, the following relation between the Presnel integrals and the tw-function 
holds [23], 



1 -f i 

C{uo) + iS{uo) = —-— ( 1 - e* 



(13) 
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Figure 1. Probability density on the t — x plane, which clearly shows the diffraction 
in time and space for an incident beam of Rubidium atoms in a fountain configuration 
[p/m = 1 cm/s, //m = 0.1 * g). The grey scale changes from light to dark as the 
function values decrease. The same color coding is used in all density plots. 



_l + i t f ft mx . 



Then, one can make use of the Faddeyeva's identity, which follows from the Cauchy 
theorem, 

w{z)+w{-z)^2e-^\ (14) 
to end up with the result 

ijj{x,t) = -e~''^'"'^'^'^w{-z), (15) 

where now 

_ 1 

^ 2 y mh^^ 2 t 
In Fig. 1 a contour diagram of the probability density on the plane t — x is plotted 
for a "fountain configuration" , namely, for the case in which the initial velocity goes in 
the direction of increasing potential. In such representation it is particularly easy to 
appreciate the diffraction in both time and space domains. 

To learn about the main features of this solution, we compare it in Fig. 2 with the 
free case. Indeed, it is amusing to rewrite (15) in terms of the Moshinsky function as 

j/^*i+i/|E ( _ ft^ p hf 

2m' h' 

which exactly reduces to the result of section 2 for / = 0. Notice that the dynamics 
is modified with respect to the free particle case by correcting both momentum and 
position with the classically expected time-dependent canonical transformation, as 
shown in Fig. 2a, up to an additional phase. In fact, the main difference is that, 
for a given u the equation of motion for the associated point of constant probability, 
and in particular for the maximum and minimum, can be deduced to be 



^|;(x, t) = e-i^+^T M a; - ^, ^, - , (17) 



pt ft^ nht 

:{t) = T; V Uo,max: (18) 

m 2m V m 



m 2m V m 
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Figure 2. a) Probability density for both free an accelerated beam at a; > with 
p/m — 5 cm/s, f /m = —0.1 cm/s^ (solid lines) and corresponding the classical density 
(dashed lines) 10 ms after removing the shutter (m = m/j^ in all the calculations), 
b) Diffraction in time for an incident beam with p/m =10 cm/s registered by a detector 
located at 5 mm from the shutter, f /m — 0.05 cm/s^. The dashed line corresponds to 
the classical density. 



where Umax = 1.2172 and Umin = 1-8725 are universal for the cut-off plane wave case 
considered here. From these equations the turning point in Fig. 2 can be inferred. 

As stated above, the telltale sign of DIT is a set of oscillations, observed after 
removing the shutter, in the probability density and which are genuinely quantum in 
nature. For diffraction in time with a linear potential, the fringe visibility 

y^ji-^ Pl^^max Pv^min (20) 

turns out to be independent of time, for the quasi-monochromatic case, as in the absence 
of the external field. 

As pointed out before the diffraction in time admits also a representation of the 
Cornu spiral. The classical probability density, given by a step function, intersects the 
quantum one at two different points. Therefore, an estimate of the width for the largest 
fringe can be carried out by considering the intersection with the classical probability 
density. The same result than without linear potential holds. 



Aa;^ 0.85a/— (21) 
V m 

[1, 2] provided that the probability is exclusively w-dependent, as should be clear from 
(11), Pmax = 1-370; Pmin = 0.778. An alternative representation which explicitly 
exhibits the diffraction pattern in time domain is shown in Fig. 2b, where the position of 
the detector is fixed and the signal of the incident beam is recorded as a function of time. 
One should notice the similarity with the profile in coordinate representation obtained by 
releasing a particle from a hard- wall trap with / = 0. This configuration was identified 
as the analogue of diffraction in time from a slit in free space [15]. Nevertheless, due to 
the dispersion entailed in (21) the pulse registered in the time domain is not symmetric. 



Matter-wave diffraction in time with a linear potential 



6 




Figure 3. Wigner function for the atom fountain configuration at the classical turning 
point {po/m = 1 cm/s, f/m = 100 cm/s^, t = 10 ms). The grey triangle in the right 
lower corner is the classically unaccessible region for which W{x,p; t) = 0. 



3.1. Wigner transform 
The Wigner function [26], 

1 

yV{x,p):=— dyip{x + yyip{x-y)e^'Py/'', (22) 

is the best known quasi joint probabihty distribution for position and momentum. 

Nowadays, it is possible to study it experimentally as has been demonstrated in a series 
of works [27]. In the context of the Moshinsky shutter, its time evolution was derived 
for the free case problem in [28]. From that result at t = 0+ or the definition (22), it 
follows that for the cut-off plane wave initial condition, the Wigner function reads 

>V(..p;po.t = Q+) = l "°'-^"'^°-'''/"i e(-x). (23) 

TT Po-P 

which clearly assumes negative values. For potentials of degree n = 2 at most, the time 
evolved Wigner function follows the classical trajectories, 

yV{x,p;t) = JJ dxodpo6[x - Xci{xo,po,t)]5[p - Pci{xo,po,t)]W{x,p;t = 0). (24) 

Alternatively one can cope with the more cumbersome calculation starting from the 
definition as in [28]. For (15) the result is 



^-{l(S + £-^)(po-^-P)}^/pt ft^ 



W{x,p;po,t) = -f ^ + ^ - X . (25) 

Fig. 3 shows the Wigner function calculated for an incident beam exactly at the classical 
turning point, being thus centered around p — 0. Note that thanks to the step function, 
it is limited to the classically accessible region of the phase space. Indeed, following [28] 
we can ask for the classical limit taking ^ ^ to find 



yVci{x,p;Po,t) ^S{po- ft-p)e ~ ■ 



(26) 
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As expected, the correct classical equation of motion x = Xci(t),p = pri{t) is recovered 
and no hint about the quantum oscillations -hallmark feature of DIT- remains. 



4. The hcird-wall trap 

In this section we consider the dynamics of a particle initially confined in a box. At 
time zero its walls are suddenly removed and a linear potential is switched on. The 
paradigmatic particle in a box (PIAB) has been recently implemented in an all optical 
fashion by Raizen et al. [13] with a Bose-Einstein condensate and experiments are in 
view with Fock states of low number of particles, N < 10 [20]. From the theoretical 
point of view, the free evolution of the wavefunction which results after shutting off the 
walls was discussed long time ago within the context of ultracold neutrons interferometry 
[14]. Later, it was reformulated by Godoy [15], who pointed out the analogy with the 
Fraunhoffer diffraction in the case of small box (compared to the De Broglie length), 
and Fresnel diffraction, for larger confinements. Quite recently , a generalization to an 
interacting many-body system was investigated in the regime of a Tonks- Girardeau gas 
[19]. 

If X[o,L]{z) is the characteristic function in the interval [0,L], let us now consider 
the initial condition given by an arbitrary excited state \ipn), 



'4'n{x,t) = W-sin (^-^j X[o,L](^)- (27) 

Using the decomposition of the characteristic function X[o,l]{x) — Q{L — x) — Q{—x), 
the ipn{x, t) term can be worked out to be 



tjjnix,t) 



2h ^'24ft 



2mfit 

L 



X 2_^a 

where we have defined 



dx e' 2w +*-2fi- _ / dx e" ^ 

-oo J —oo 



(28) 



ahmr ft , , 

Carrying similar steps to the ones described for the previous section, one finally ends 
up with 



ay L 

where 



a=± 



h 2ftt w{-UaL) - 2« w{-UaG) 



(30) 



1 + i I t ( m(x - L)\ 1 + i I t ( mx\ 

= ^ l^" — r- j ' = ^ - ~) -^^'^ 

In an atomic fountain, trapped atoms arc launched vertically with the help of the optical 
molasses technique. PIAB eigcnstatcs are a good approximation to the actual state of 
a box subjected to a linear potential, whenever the box is small and the particle light. 
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t{ms.) 



Figure 4. Time evolution of the density profile for an eigenstate of the hard-wall 
trap, with n = 10, q = cm/s, L = 80 fim and f/m = 0.49 cm/s^. Bifurcation for 
t > tn=io — 137 ms takes place in analogy with free case. The expectation value {x{t)) 
reproduces the classical trajectory. 



Indeed, within first order perturbation theory the coefficients Cnk responsible for the 
corrections on the n-th eigenstate (|0n'*) = + ^kj^n^rik\4>^k^)) due to the linear 

potential are given by 

^nk - . (32) 

If the perturbation of the initial state by the linear potential is significant, it can 
always be written as a linear combination of the eigenstates of the free Hamiltonian 
within the box, namely, {(pni^) — ^^2/ L sm{n7rx / L)\n e N}. We must therefore look 
for the evolution of a PIAB state released at time t — and launched with a given 
momentum q. Explicitly, it is described by (30), with the redefinition 
_ ahmr ft 

Pa = q + — y • (33) 

Equation (30) can model too the output coupling mechanism of an atom laser. In such 
devices, an atom is initially in the lasing mode, confined by the cavity mirrors which can 
be simulated by two infinite wall (reducing the problem to that of a PIAB). Through a 
Raman transition, this atom evolve to a nontrapped state, which amounts to shutting 
off the walls. During the transition and due to the emission of a photon the atom is 
kicked with a given momentum, namely, q [29] . 

Figure 4 shows the time evolution of the density profile for the tenth eigenstate of 
the initial hard-wall trap. A most relevant feature is that at variance with the well-known 
harmonic trap case, the profile bifurcates into two main branches after the semiclassical 
time tn = mL/{2pn) = mL"^ / {2mili) . This behaviour has been recently reported for the 
free particle case by the authors [19] and holds for any excited state for both even and 
odd quantum number n > 1. The cloud falls under acceleration following ((ilj\x{t)\ip) the 
classical trajectory, even though for any t > tn the probability is vanishing along it. 
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5. Tonks gas dynamics in a linear potential 

Under strong radial confinement and at low densities and temperatures an ultracold 
atomic gas becomes effectively one dimensional. In particular, whenever the thermal 
and zero-point energies are lower than the transversal excitation quantum, the Tonks- 
Girardeau (TG) regime is reached. The effective interactions are then so strong 
that the gas behaves like a system of hard-core impenetrable bosons. Such regime 
has been obtained in several experiments [30] and a full quantum description is 
possible thanks to the Fermi-Bose mapping theorem [31]. According to it, the many- 
body wavefunction of a TG gas can be obtained from the one of a zero-spin free 
fermionic gas, by applying the "antisymmetric unit function" A = ni<j<A:<Af ^9''^i^k — 
Xj), as iJjb{xi, . . . ,XN;t) = A{xi, . . . ,xj^)ipp{xi, . . . ,XN;t). Indeed one can state 
\iPb{xi, . . . jXN'jt)]"^ = \iPf{xi, . . . ,XN',t)\'^, for the mapping is involutive. The dual 
system to the TG gas, being the ideal Fermi system, is built up as a Slater determinant 



Moreover, thanks to the fact that under unitary time evolution, orthonormality between 
the initial eigenstates of the trap is preserved, it follows that any local correlation 
function can be obtained in close-form. For example, the calculation of the density 
profile is greatly simplified to 



The free expansion dynamics of a TG gas from a hard- wall trap was recently considered 
in [19], where a dynamics much more complicated than for the harmonic confinement 
was observed. Taking advantage of the results obtained in the previous section wc next 
generalize the free expansion results and deal with the following experiment. Consider 
a TG gas initially confined in a hard wall trap. At zero time it receives a momentum 
kick, momentum at which the shutter is off and the linear potential rumped up. The 
density profile is found combining (30) and (35). 

Figure 5 shows the density profile of the expanding cloud in space and time falling 
under the action of a linear potential. For short times t < t^^ there is a well-defined 
pattern where the number of maxima equals that of particles. After a transient regime 
where the visibility of the peaks varies in a non-trivial way, it is gradually lost for f > tN- 
In the later regime, the ballistic expansion is established in agreement with the force-free 
case [19]. 

6. Conclusions 

Diffraction in time has been generalized in the presence of a linear potential. Moreover, 
the dynamics of particles released from a hard-wall trap under a constant force has 
been studied, showing a bifurcation after a time tn — mL^ /{2n7rh), in analogy with 




(34) 



g{x, t) = N \'iIjb{x,X2,. . . , xn] t)\'^dx2 ■ ■ ■ dxn = ^ |0n(a;, t) |^. (35) 



n=l 
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Figure 5. Time evolution of the density profile for Tonks- Girardeau gas , with N = 10, 
q = cm/s, L = 80 yum and f/m = 0.49 cm/s^. The interference pattern is washed 
out for t > tjv=io — 137 ms. 

the free case [19]. This is a novel feature exclusively associated with hard- wall traps. 
Note that for the case of harmonic confinement, the evolution follows a scaling law of 
coordinates, lacking therefore any transient effect [33]. In addition, the expansion of 
a Tonks- Girardeau gas has been exactly solved using the Bose-Fermi map. The cloud 
is shown to exhibit an interference pattern which is lost for t > t^. Finally, we point 
out that any time dependence on the linear potential can be easily included using the 
general map obtained in [34] . In such a way, one can account for a smooth ramping up 
of the linear potential. 
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